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ABSTRACT 

The a Cen binary system is a well-known stellar system with very accurate obser- 
vational constraints to structure of its component stars. In addition to the classical 
non-seismic constraints, there are also seismic constraints for the interior models of a 
Cen A and B. These two types of constraint give very different values for the age of 
the system. While we obtain 8.9 Gyr for the age of the system from the non-seismic 
constraints, the seismic constraints imply that the age is about 5.6-5.9 Gyr. There 
may be observational or theoretical reasons for this discrepancy, which can be found 
by careful consideration of similar stars. The a Cen binary system, with its solar type 
components, is also suitable for testing the stellar mass dependence of the mixing- 
length parameter for convection derived from the binaries of Hyades. The values of 
the mixing-length parameter for a Cen A and B are 2.10 and 1.90 for the non-seismic 
constraints. If we prioritize to the seismic constraints, we obtain 1.64 and 1.91 for 
a Cen A and B, respectively. By taking into account of these two contrasting cases 
for stellar mass dependence of the mixing-length parameter, we derive two expres- 
sions for its time dependence, which are also compatible with the mass dependence 
of the mixing-length parameter derived from the Hyades stars. For assessment, these 
expressions should be tested in other stellar systems and clusters. 

Key words: stars: interior - stars: evolution - stars: individual: a Cen ~ stars: 
abundances - binaries: visual - stars: oscillations 



1 INTRODUCTION 

The observable quantities of a star, such as luminosity, ra- 
dius and effective temperature, mostly depend on its mass. 
Chemical composition, rotation and age, however, have a 
second role in these quantities. The time dependence of the 
observable quantities are very weak. Therefore, no more or 
less precise age for stars available, particularly for late-type 
stars, and this is an essential quantity for our understanding 
of the evolution of both the far and near Universe, However, 
it is possible to derive the ages of stars from their seismic 
properties (Christensen-Dalsgaard, 1988). In this respect, 
the solar-like components of a Centauri (a Cen) are excel- 
lent targets, as their seismic properties are determined by 
ground-based spectrographic observations. The aim of the 
present paper is to test the theory of stellar evolution by 
constructing models for the internal structure of these stars, 
using all the available constraints. 

The internal structures of ct Cen A and B have been 
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considered theoretically in many papers. Recent studies 
by Miglio & Montalban (2005), Eggenberger et al. (2004), 
Thoul et al. (2003), Morel et al. (2000) and Guenther & De- 
marque (2000) have analysed the structure and evolution of 
these stars in detail. The detection of their seismic proper- 
ties have increased the scientific interest in these stars. The 
large and small frequency separations of a Cen A and B 
have been derived observationally from their p-mode oscil- 
lations by Bouchy & Carrier (2002), and Carrier & Bourban 
(2003), respectively, using spectrographic methods. In ad- 
dition to the constraints from the fundamental properties 
of the stars, the interior models of these stars should also 
satisfy these seismic constraints. Recently, Kjeldsen et al. 
(2005) also observed a Cen B, and measured the frequen- 
cies, which turned out to be quite different from the previous 
results of Carrier & Bourban (2003). Therefore, the system 
deserves to be reconsidered by taking this new measurement 
into account. 

The masses (Pourbaix et al. 2002) and radii (Kervella 
et al. 2003) of a Cen are well known from the cited observa- 
tions. Although the a Cen system has already been widely 
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investigated, no agreement lias yet boon reached about the 
precise value of the effective temperatures (Teff) and ele- 
ment abundances in the atmospheres of these stars. These 
quantities are determined by spectroscopic methods. How- 
ever, using the photometric observation of Bcsscl (1990) as 
an alternative, which has unfortunately been neglected in 
studies on a Cen, we have alternatively derived the effective 
temperatures and the heavy element abundances of both 
components. 

The effective temperature of a model is a function of the 
mixing-length parameter for convection. Yildiz et al. (2006) 
have confirmed a very definite relationship between the 
mixing-length parameter and the stellar mass by fitting the 
models of late type components of some binaries in Hyades 
to the observations: a = 9.19(M/Mq - 0.74:)° °^^ - 6.65. 
However, controversial results are obtained from studies on 
binaries by different investigators: while Lastennet et al. 
(2003), Ludwig & Salaris (1999), and Lebreton et al. (2001) 
have stated that a is an increasing function of the stellar 
mass in their studies on some other binaries, Miglio & Mon- 
talban (2005), Eggenberger et al. (2004), Morel et al. (2000) 
have confirmed that ob is greater than oa for the a Cen 
binary system. In these studies and others, at least for the 
favoured models, however, we note that the model luminos- 
ity is less than the observed value for a Cen A, and that 
the opposite is true for a Cen B. This point needs to be 
explained. However, the results of the hydrodynamic sim- 
ulations of convection arc in agreement with the decrease 
a. Thus, there may be no single relationship between the 
mixing-length parameter and the stellar mass and we should 
therefore search for a more general expression for a than 
that given above. 

According to the above-mentioned studies on the inter- 
nal structure of a Cen A and B, the age of the system is 
between 4.85-7.6 Gyr, at least for the favoured or typical 
models. The lower limit is given by Thevcnin ct al. (2002) 
and the upper limit is from Gucnthcr & Dcmarque (2000). 
Miglio & Montalban (2005), however, find the age of the 
system to be 8.9 Gyr from the models for the non-seismic 
constraints (see Section 4.1). In other words, the seismic and 
the non-seismic constraints lead to very different ages for a 
Cen. 

The remainder of this paper is organized as follows. In 
Section 2, we summarize the observed chemical composition, 
the seismic and the fundamental properties of a Cen A and 
B. The description of models are given in Section 3, and 
the results of the model computations are presented and 
discussed in Section 4. Finally, we give concluding remarks 
in Section 5. 



2 OBSERVED PROPERTIES OF a CEN 

2.1 Chemical abundance of a Cen A and B 

Many studies have been done on the abundance determi- 
nation of chemical species of a Cen from the spectra of its 
components. Most of these studies deal with the abundance 
of Fc and/or similar heavy elements, which arc not abundant 
in stars. The most complete studies-in the sense that they 
also deal with the most abundant heavy elements in normal 
stars such as oxygen, carbon and nitrogen-have been cwcar- 



ried out by Neuforgc-Vcrheecke & Magain (1997) and Feltz- 
ing & Gilmore (2001, hereafter FG2001). The abundances 
of 17 elements have been obtained by Neuforge-Verheecke & 
Magain (1997). They find from the spectra of a Cen A that 
overabundance relative to solar is 0.25 dcx for Fe and 0.21 
for the most abundant heavy element oxygen, and a general 
overabundance of 0.24 dex. They also confirm that there is 
no significant difference between the chemical compositions 
of a Cen A and B, but a Cen B is likely more metal rich 
than a Cen A. Using their equivalent widths, FG2001 have 
obtained very similar results for the abundance of 13 heavy 
elements, including oxygen. 

By considering Mg, Ca, Si, Ti, Cr and Fc elements, 
Doyle ct al. (2004) have found that the mean overabundance 
of q: Cen A is 0.12±0.06 dex. Contrary to this result, Ecuvil- 
lon et al. (2006) computed the average value of the oxygen 
abundance of a Cen A and B as [O/H]=-0.12 and -0.06, 
respectively, from triplet lines by applying non-local ther- 
modynamical equilibrium effects. Such contradictions in the 
abundance determination a Cen A and B, and the recent 
revisions on the solar compositions (Asplund et al. 2005) 
make us cautious about seeing the abundance determina- 
tion of stars as an issue that has been conclusively studied. 

The mixture in FG2001 is taken as the chemical compo- 
sition of the OPAL opacity tables (Iglesias & Rogers 1996) 
used in the construction of models for a Cen A and B. For 
the chemical species that are not observed by FG2001, such 
as Nc, wc adopt the value of 0.25 dcx for their abundance 
relative to the Sun. For comparison, we also construct mod- 
els with the recent solar mixture given by Asplund et al. 
(2005, hereafter AGS2005). 

2.2 Fundamental properties of a Cen A and B 

For our better understanding of stellar structure and evolu- 
tion, a very accurate knowledge of stellar masses is required. 
The determination of masses of a Cen A and B from ob- 
servations has been the subject of many papers. The most 
recent determination is that of Pourbaix et al. (2002). Ac- 
cording to their findings, the masses of a Cen A and B are 
as Ma = 1.105 ± 0.007 M® and Mb = 0.934 ± 0.007 Mq, 
respectively. The radii of the components have been mea- 
sured by Kervella et al. (2003) as Ra = 1.224 ± 0.003 R© 
and Rb = 0.863 ± 0.005 R©. The uncertainties in radii are 
computed by taking into account only the uncertainty in the 
parallax. Many investigators have determined the effective 
temperatures from the spectra of a Cen A and B is done 
by many investigators (see, for example. Morel et al. 2000 
and references therein). According to the results of studies 
based on spectra, while the effective temperature of a Cen 
A is between 5830 K and 5720 K, the range of effective 
temperatures of a Cen B is 5250-5325 K. 

In the literature, apparent magnitude (V) and colour 
(B-V) of a Cen A and B are, in general, taken from The 
Bright Star Catalogue (Hoffleit & Jaschek 1982): Va = 
-0.01, Vb = 1.33, (B - V)a = 0.71 and (B - V)b = 0.88. 
Adopting the parallax given in Soderhjelm (1999), tt = 
747.1 ± 1.2 mas, we find the absolute magnitudes of a Cen 
A and B as 4.358 and 5.697, respectively. Using the tables 
of Lejeune et al. (1998, hereafter LCB1998) for colours and 
bolometric correction, prepared for different heavy element 
abundances, we obtain the required luminosities of a Cen 
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Table 1. The fundamental properties of a Cen A and B. Using visual magnitudes (HofHeit &c Jaschek 1982; Bessel 1990) and the 
observed masses and radii of a Cen A and B, we find bolometric corrections from the tables of LCB1998 and Besl998, and then their 
luminosity for selected values of Z. The results given in the last two rows are obtained by also fitting the colours of a Cen A and B. In 
this way, we also determine the current heavy element abundance of each star. 



Star 


La/I/0 




MboI 


BC 


U -B 


B-V 


My 


[Z/Ze] 


Table(BC) 


a 


Cen A 


1.552 


5831 


4.2728 


-0.0852 


0.2042 


0.6583 


4.3579 


0.25dex 


LCS1998 


a 


Cen A 


1.555 


5834 


4.2707 


-0.0876 


0.1837 


0.6509 


4.3583 


0.20dex 


LCB1998 


a 


Cen A 


1.585 


5862 


4.2499 


-0.1081 


0.0721 


0.6077 


4.3580 


O.OOdex 


LCB1998 


a 


Cen A 


1.520 


5801 


4.2954 


-0.0623 


0.1578 


0.6642 


4.3577 


O.OOdex 


Besl998 


a 


Cen B 


0.511 


5260 


5.4796 


-0.2174 


0.5922 


0.8559 


5.6969 


0.2bdex 


LCS1998 


a 


Cen B 


0.511 


5261 


5.4789 


-0.2176 


0.5690 


0.8486 


5.6966 


0.20dex 


LCB1998 


a 


Cen B 


0.515 


5271 


5.4705 


-0.2261 


0.4258 


0.8037 


5.6966 


O.OOdex 


LCB1998 


a 


Cen B 


0.497 


5225 


5.5082 


-0.1885 


0.4746 


0.8552 


5.6968 


O.OOdea; 


Sesl998 


a 


Cen A 


1.544 


5824 


4.2784 


-0.1019 


0.1276 


0.6333 


4.3803 


Q.midex 


LCB1998 


a 


Cen B 


0.507 


5250 


5.4875 


-0.2227 


0.5349 


0.8398 


5.7102 


0.l2bdex 


LCB1998 



A and B for selected values of Z. In Table 1 the first three 
rows for a Cen A and B are for the heavy element abun- 
dance of 0.25 dex and 0.20 dex and for the solar metallicity, 
respectively. For each metal abundance, while wc obtain rea- 
sonable values for the absolute magnitudes from the tables 
in comparison with the observed values given above, for the 
colours, there is no agreement between theoretical and ob- 
servational results. In the fourth and the eighth rows, we 
derive the fundamental properties of a Cen A and B, re- 
spectively, using the tables of Bessel, Castelli & Plez (1998, 
hereafter Besl998) for colour and bolometric correction, pre- 
pared with the solar composition. While the luminosities 
of a Cen A and B found from Besl998 are smaller than 
the corresponding luminosities computed from the tables of 
LCB1998, the colours computed from the two tables axe sig- 
nificantly less than the observed colours for all the selected 
metallicities. 

In Bessel (1990), however, the apparent magnitudes and 
the colours of (B-V) a Cen A and B are given as Va = 0.01, 
Vb = 1.34, {B - V)a = 0.633 and {B - V)b = 0.84. 
These apparent magnitudes give absolute magnitudes of 
Mv,A = 4.377 and Mv,b = 5.707. Using the above mass, 
radius, Mv and B — V each component we derive the lu- 
minosity, effective temperature and also the current value of 
the surface heavy element abundance (Zs) from the tables of 
LCB1998. Their tables are prepared with different values of 
Z, which enable us to determine Zs - By changing Z and the 
luminosity of the component stars, we obtain the observed 
absolute magnitude and (B-V) of each star. The fundamen- 
tal properties of a Cen A and B we have derived are listed 
in the last two rows of Table 1, respectively. The current 
values of the surface heavy element abundances of a Cen A 
and B relative to the Sun are found as 0.074 dex and 0.125 
dex, respectively. The value of Z depends on what value is 
adopted for Zq in the tables of LCB1998. For Zq = 0.02, 
Zs,A = 0.0237, Zs,B = 0.0267. While the luminosities we 
find are slightly greater than those given by Eggenberger 
et al. (2004), the effective temperatures of a Cen A and B 
are found as 5824 K and 5250 K, respectively, and are in 
very good agreement with the values found from spectra by 
Neuforge-Verhccckc & Magain (1997). 

Indeed, the measurements given in Bessel (1990) and 
HofHeit & Jaschek (1982) are very close to each other; the 
remarkable differences that exist between the colours of 
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Figure 1. Observational seismic HR diagram for the Sun (0), a 
Cen A (filled circle) and B (star), r; Boo (box), f3 Vir (diamond) 
and Procyon A (circles). The seismic data are taken from Chaplin 
et al. (1999), Bouchy & Carrier (2002), Kjcldscn ct al. (2005), 
Carrier et al. (2005a, b) and Martic et al. (2004). Do is computed 
from the small separations 51^02 between the observed frequencies. 
For the Sun and a Cen B, 5:^13 is also available and is thus also 
used in the computation of Dq of these stars. 



each star and the measurements of Bessel (1990) are com- 
patible with the theoretical results. Therefore, we adopt the 
photometric constraints in Bessel (1990) in the calibration 
process of models for a Cen A and B. 



2.3 Seismic properties of a Cen A and B and an 
age estimation for the system 

The most important constraints on the stellar evolution the- 
ory from the observations of solar-like oscillations are small 
{Svni) and large (Az/„j) separations between the observed 
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frequencies. The small separation between the frequencies 
(Christensen-Dalsgaard 1988) is defined as 

5l>nl = ^71,1 — ^'n-l,i+2- (1) 

In place of this, Do is frequently used as an average effect of 
oscillations with different values of degree I. Wc compute Do 
from 5vnO'- Do = Svno/Q (Kjeldsen et al. 2005). If frequencies 
for / = 3 modes are observed, Do is computed from SvnO and 
SVr^V.Do = (6uno/6 + Siyni/10)/2. 

The large separation between the frequencies is given 
below: 

Ai/„i = Un,i — iy„-i,i. (2) 

While the small separation is a sensitive function of physi- 
cal conditions in central regions where the rmclcar evolution 
proceeds, the large separation is a measure of the mean den- 
sity. 

The large and small frequency separations of a Cen A 
are determined by Bouchy & Carrier (2002) as 105.5 ± 0.1 
fiRz and 5.6 ± 0.7 ^iHz (Do = 0.93±0.12 /xHz), respectively. 
Carrier & Bourban (2003) find that SfnO = 8.7 ± 0.8 /xHz 
and Avno = 161.5±0.1 fjHz for a Cen B. Recently, Kjeldsen 
et al. (2005) also obtained the seismic properties of a Cen B 
using the two-site observations. According to their findings, 
while the large separation between the frequencies is very 
similar to that of Carrier & Bourban (2003), the small sep- 
aration between the frequencies is Sfno = 10.14 ± 0.62 /xHz 
and significantly larger than the previously found separa- 
tion (8.7 ±0.8 /iHz) by Carrier & Bourban (2003). Kjeldsen 
et al. (2005) also observed the oscillations with / = 3 and 
the small sepaxation between the frequencies of these oscil- 
lations and oscillations with Z = 1 is Sfni = 16.73 d= 0.65 
/iHz. From Svni.) and Si^ni, Do — 1.68 ± 0.08 /iHz is found. 

Fig. 1 shows and observational seismic HR diagram for 
several stars for which Dq and Av^i are available. While the 
filled circle and star represent a Cen A and B, respectively, 
© shows the position of the Sun. Do and Avni of the Sun 
arc computed from the data of Bison group (Chaplin et al. 
1999). The seismic data of rj Boo (box), /3 Vir (diamond) 
and Procyon A (circle) are taken from Carrier, Eggenberger 
& Bouchy (2005), Carrier et al. (2005), and Martic et al. 
(2004), respectively. The Do values of a Cen B and the Sun 
arc very close to each other. This means that the central hy- 
drogen abundance of a Cen B {Xc = 0.34) is nearly equal to 
the solar value (see figure 1 in Christensen-Dalsgaard 1988). 
Thus, the age of a Cen B is about half of its nuclear time 
scale (tnuc). From this confirmation, we can estimate the 
age of a Cen using a simple approach. Because the nuclear 
time-scale is inversely proportional to M^'^ (Bohm-Vitense 
1992), 

(j\/r \ 2.5 
^ j tQ = 5.5 Gyr. (3) 

The basic observational properties of a Cen A and B 
discussed in this section are given in the last two rows of 
Table 2. 



3 DESCRIPTION OF THE MODELS 

The characteristics of our code have already been described 
by Yildiz (2000, 2003, see also references therein), and there- 
fore we shall not go into details. Our equation of state uses 




0.2 I ' ' ' ' ' 1 

6000 5800 5600 5400 5200 5000 4800 

Teff(K) 

Figure 2. Evolutionary tracks of models (NOS model) for a Cen 
A (thick line) and B (thin line) in the HR diagram. While the 
diamond represents the data wc derived from photometric data 
in Bessel (1990), the circle shows the data given by Eggenberger 
et al. (2004) 

the approach of MihalEis et al. (1990) in the computation 
of the partition functions. The radiative opacity is derived 

from Iglesias & Rogers (1996), completed by the low tem- 
perature tables of Alexander & Ferguson (1994). For the 
nuclear reaction rates, we use the analytic expressions given 
by Caughlan & Fowler (1988), and we employ the standard 
mixing-length theory for convection (Bohm-Vitense 1958). 
For the diffusion of chemical species, the routines of Thoul, 
BahcaU & Loeb (1994) are used. 



4 RESULTS AND DISCUSSIONS 

4.1 Models satisfying the constraints on the 
fundamental properties of a Cen A and B 

In order to find the five unknowns (i.e. X, Z, t, oa and as), 
we write down four equations for La, Lb, Ra and Rb- The 
fifth equation for a solution is for the surface value of the 
(Z/X) ratio for a Cen A. For {Z/X)a,s = 0.0237, the values 
of five unknowns and the models of a Cen A and B using 
these values are given in the first two rows of Table 2. The 
heavy element mixture of the models is taken from FG2001. 
According to the solution for non-seismic constraints (NOS 
model), the age of the system is 8.88 Gyr and a Cen A 
is very close to the end of its main-sequence (MS) lifetime. 
Because the seismic constraints are not used in this solution, 
the seismic properties of models are not in agreement with 
the observations: small separations between the frequencies 
of the models are significantly smaller than the observed 
separations for both stars. In the last two rows of Table 2, 
the observational results are listed for comparison. In Fig. 
2 the models are plotted in a theoretical HR diagram. The 
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Table 2. Model properties of a Cen A and B. <5j^o2, Suis, Do and Ai/q are In units //Hz. The uncertainty of 61/02 and 61/13 Is about 0.7 
//Hz. The error In Auo Is 0.1 //Hz. The uncertainty of Do Is 0.12 //Hz for a Cen A and 0.08 //Hz for a Cen B. 



Star 


L/Lq 


R/Rq 


Teff 


^0 


Zo 




{Z/X)s 


a 


t(109y) 


61/02 


6vi3 


Do 




cc 


MODEL 


A 
B 


1.545 
0.507 


1.224 
0.864 


5822 
5246 


0.703 
0.703 


0.0328 
0.0328 


0.0237 
0.0267 


0.0316 
0.0364 


2.10 
1.90 


8.88 
8.88 


2.7 
6.6 


15.8 


0.45 
1.30 


106.5 
164.6 


r GrzOOl 
FG2001 


JMUo 
NOS 


A 

B 


1.545 

0.584 


1.224 

0.863 


5822 

5436 


0.669 

0.669 


0.0322 

0.0322 


0.0237 

0.0271 


0.0316 

0.0388 


1.64 

1.91 


5.70 

5.70 


5.57 

9.15 


18.6 


0.93 

1.68 


107.4 

164.6 


FG2001 

FG2001 


SIS 

SIS 


B 


0.506 


0.862 


5248 


0.669 


0.0322 


0.0271 


0.0389 


1.58 


5.70 


9.50 


19.0 


1.74 


163.5 


FG2001 


5/5914 


A 
B 


1.543 
0.584 


1.223 
0.863 


5822 
5441 


0.676 
0.676 


0.03 
0.03 


0.0198 
0.0271 


0.0316 
0.0388 


1.65 
1.93 


5.57 
5.57 


5.64 
9.29 


18.4 


0.94 
1.70 


106.8 
164.6 


FG2001 
FG2001 


5/53 
5/53 


A 
B 


1.549 
0.569 


1.223 
0.863 


5828 
5400 


0.713 
0.713 


0.023 
0.023 


0.0161 
0.0192 


0.0211 
0.0259 


1.71 
1.98 


5.90 
5.90 


5.49 
9.31 


18.5 


0.91 
1.70 


106.7 
164.4 


FG2001 
FG2001 


5/523 
5/523 


A 
B 


1.544 
0.572 


1.223 

0.863 


5828 

5409 


0.697 
0.697 


0.0233 
0.0233 


0.0159 
0.0192 


0.0213 
0.0265 


1.67 
1.96 


5.70 
5.70 


5.61 
9.38 


18.6 


0.93 
1.71 


106.8 
165.0 


AGS2005 
AGS2005 


SISAGS 
SISAGS 


A 
B 


1.544 
0.507 


1.224 
0.863 


5824 
5250 






0.0237" 
0.0267" 


0.0364'' 
0.0399= 






5.6 
10.14 


16.73 


0.93 
1.68 


105.5 
161.5 




obs 
obs 





"Derived from the photometric data of Bcssel (1990) using the tables for colours and bolometric corrections of LCB1998. 
) Derived from abundances found by FG2001. For the elements not observed by them, an average overabundance of 0.25 dex relative 

to solar is assumed. 

°) Derived from assuming an overabundance of 0.04 dex relative to a Cen B (observed for Al and Si, for example). 




observed positions of a Cen A and B derived from the pho- 
tometric measurements (diamond) arc also soon. For com- 
parison, the data given by Eggenberger et al. (2004) are 
also plotted (circle). The uncertainties in the observational 
luminosities and the effective temperatures are taken from 
Eggenberger et al. (2004). According to the NOS model, 
while a Cen A is very close to the end of its MS lifetime, a 
Cen B is at about half of its MS lifetime. The value of the 
mixing-length parameter found for a Cen A is greater than 
that for a Cen B. This result is in good agreement with the 
results given by Yildiz et al. (2006) . 

The ago of those models, 8.88 Gyr, also found by Miglio 
& Montalban (2005), is almost the same as the age of the 
galactic thin disk (8.8 Gyr) obtained by del Peloso et al. 
(2006). Their age determination is based on the abundance 
ratio of Th and Eu elements in 26 stars including a Cen. 



4.2 Calibration of models with seismic constraints 

As discussed in the previous subsection, the small separa- 
tions of NOS models of a Con A and B are significantly 
smaller than the observed values. Therefore, models satis- 
fying the seismic constraints must have a smaller age. How- 
ever, models of a Cen A and B require a high value for 
age, such as 8.88 Gyr, in order to fit the observed difference 
between the luminosities of a Cen A and B. Because there 
is no simultauoous solution for La, Lb and the seismic con- 
straints, wo therefore take into account the luminosity of one 
component only. If we leave aside the luminosity of a Cen B 
and use all the remaining constraints including seismic mea- 
surements, wo find a significantly smaller age for the system 
(SIS model) than predicted by NOS models. The SIS models 
of a Cen A and B are given in the third and fourth rows 
of Table 2, respectively. In contrast to the NOS models, the 
value of the mixing-loiigtli parameter found for a Cen A is 
smaller than that of a Cen B. However, the reason of this 
result could be that the luminosity of the SIS model for a 
Cen B is 15 per cent greater than the observed luminosity. 
If this excess is a result of observational uncertainty 
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Figure 3. Theoretical seismic HR diagram for a Cen A and B 
(SIS model). While tlie dotted lines show the evolutionary tracks 
of the stars, the solid lines show the isochrones. The solid line at 
the top corresponds nearly the ZAMS (O.lGyr). The solid lines 
below it show 1, 2, 3, 4, 5 and 5.7 Gyr. For comparison, the 
observed positions of the Sun(0) and a Cen A (filled circle) and 
B (diamond) are also plotted. 

about the mass of a Cen B, then its required mass is 0.914 
Mq. The properties of the model with this mass are given 
in the fifth row of Table 2. The mixing-length parameter of 
this model is less than the corresponding model of a Cen A. 

Because the abundance determination of stars is not an 
issue that has been pursued exhaustively, we also construct 
models for a Cen A and B with two different Z. For Z = 0.03 
(sixth and seventh rows in Table 2) and Z = 0.023 (eighth 
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and ninth rows), putting aside the differences between the 
initial hydrogen abundances, we find very similar models. In 
order to see the effect of the heavy element mixture, we also 
obtain models with the recent solar heavy element mixture 
(AGS2005) in place of that of FG2001. 

In Fig. 3, the time run of the models of a Cen A and B 
(SIS models) are plotted in a seismic HR diagram. The filled 
circle and diamond show a Cen A and B, respectively, and 
the Sun is represented by ©. While the dotted lines mark the 
evolutionary tracks, the solid lines show the isochronos. The 
isochrone line at the top is for 0.1 Gyr and the next line is for 
1 Gyr. The increment between the other isochrone lines is 1 
Gyr, except the last one. The isochrone line at the bottom 
is for 5.7 Gyr. 



4.3 On the alternative seismic HR diagrams 

Seismic and non-seismic constraints impose very different 

ages on the models of a Con A and B. For the small separa- 
tion between the frequencies, however, oxburgh & Vorontsov 
(2003) and Mazumdar (2005) propose alternative expres- 
sions for a better representation of physical conditions in 
central regions of the late-type stars. Roxburgh & Vorontsov 
(2003), for example, propose to use 

doi(n) = \ {Vn-lfi - 4t'n-l,l + QVnfi " 4fc'„,l + Z/„+l,o) (4) 
o 

in place of Do (or 5i'o2)- They also confirm that the ratio of 

small to large separations (roi = doi(n)i/Ai'„;) for acoustic 
oscillations is more sensitive to physical conditions in the 
stellar core. From the observed seismic data of a Cen A 
and B, we find that roi is 0.018 and 0.025, respectively. The 
uncertainty in doi(n) is, in principle, the same as for the 
small separations. Then, the uncertainty in roi is 0.007 for 
a Cen A and 0.004 for a Cen B. For the SIS model, we 
compute the ratio roi = do\{n)i / Avni'- roi = 0.045 for a 
Cen A and roi = 0.012 for a Cen B. Very similar results are 
also obtained for the other models listed in Table 2. While 
the value of roi for the model of a Con A is significantly 
larger than the observed roi and does not show an explicit 
time dependence, the observed roi of a Cen B is close to 
the zero-age main-sequence (ZAMS) value of the SIS model 
rather than the value at 5.7 Gyr. 



4.4 Variability of the mixing-length parameter 
with time 

As stated above, Yildiz et al. (2006) have confirmed a re- 
lationship between the mixing-length parameter for convec- 
tion and the stellar mass from the late-typo stars of the 
Hyades open cluster. According to this relationship, for a 
Cen A and B, qa = 2.06 and ae = 1.78. A similar relation- 
ship is found for the models of a Con A and B fitted into 
the non-seismic constraints of the system. The mixing-length 
parameters for these models are oa = 2.10 and as = 1.90. 
These parameters are qualitatively in very good agreement 
with the results of Yildiz et al. (2006) from the binaries of 
the Hyades open cluster: the higher the mass, the greater 
the mixing-length parameter (case 1). 

However, this relationship may have different forms as 
the stars evolve away from the ZAMS (see below). For all 




0.7 0.8 0.9 1.0 1.1 

M/Mq 



Figure 4. The solid line represents the stellcir mass dependence 
of the mixing-length parameter given by Yildiz et al. (2006) for 
the Hyades stars. W'e try to derive more general expressions for a 
in terms of the structure of the convcctivc zones rather than the 
stellar mass. For two typical cases, we plot /i and /2 as functions 
of some physical quantities of the convective zone (such as den- 
sity and temperature at the base of convective zone and at the 
surface). See the text for the definitions of /i and /2. 
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Figure 5. The time dependence of the mixing-length parameter 
for a Cen A predicted by /i (circles) and /2 (stars). 
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Figure 6. The time dependence of the mixing- length parameter 
for a Cen B predicted by fi (circles) and /2 (stars). 



models satisfying the seismic constraints, for example, con- 
trary to case 1, QfA is less than qb (case 2). In case 2, at least 
for a Cen A, a may be a decreasing function of time. Nei- 
ther case 1 nor case 2 gives the same values for qa and as- 
This means that a is not a simple function of stellar mass 
and the mass dependence found by Yildiz et al. (2006) (see 
their figure 4) is a special case of a general situation about 
the ZAMS. So, we describe a as a function of some quanti- 
ties pertaining to the convective zone of each model rather 
than stellar mass. In Fig. 4, we plot two such functions, /i 
(circles) and /2 (stars), with respect to stellar mass: 



/i = 2.5 - pbc 



2.7 

fe.bcz 



1.4 

fs.ph 



_ 3^ 

„0.8 



V 3.8 ; 



(5) 



(6) 



Here, the subscripts bcz and ph mean that the correspond- 
ing quantity pertain to the base of convective zone and the 
photosphere, respectively, rbcz is the radius of the base of 
the convective zone and Afbcz is the mass inside the sphere 
with this radius, rbcz and Afbcz arc in solar units. The 
solid line in Fig. 4 represents the fitting curve a{M) = 
9.19(M/M0 -0.74)" °5^ -6.65 derived by Yildiz et al. (2006). 
The three functions predict very similar values of a for a 
given mass. However, there are very important differences 
between these three functions. While a{M) is constant in 
time, /i and /2 are implicitly functions of time because the 
structure of the convective zone changes in time. 

In Fig. 5, /i (circles) and /2 (stars) are plotted as a 
function of time for a Cen A (SISZ23 model). While /i is 
an increasing function of time, /2 is a decreasing function. 
We also see how /i and /2 vary in time for a Cen B. In Fig. 
6, we see that both /i (circles) and /2 (stars) are increasing 



functions of time and predict almost the same value for a 
Cen B for a given time. 



5 CONCLUSIONS 

The masses and radii of the components of a Cen, which 
arc essential quantities for stellar modelling, are accurately 
determined from the observations. Using the photometric 
measurements of Bessel (1990) and the tables for colours 
and bolometric corrections of LCB1998, we find the abso- 
lute magnitudes of a Cen A and B and then their luminosi- 
ties and surface heavy element abundances: La = 1.544Lq, 
Lb = O.5O7L0, 0.074 dex and 0.125 dex for the metallic- 
ities of a Cen A and B, respectively ( Za = 0.0237 and 
Zb = 0.0267, taking the metallicity of the tables of colours 
and bolometric correction as Z=0.02). This difference be- 
tween the metal abundances of a Cen A and B can be inter- 
preted as the result of the diffusion process (settling), which 
is faster at the bottom of convective zone of a Cen A than 
a Cen B. 

The calibration of models of a Cen A and B into these 
uon-soismic constraints yields the age of the system as 8.88 
Gyr. The mixing-length parameters for these models are 
a A = 2.10 and as = 1.90. These values of a are qualita- 
tively in very good agreement with the results of Yddiz et 
al. (2006) from the binaries of the Hyades open cluster: The 
higher the mass, the greater the mixing-length parameter. 
However, this relationship may have different forms as the 
stars evolve away from the ZAMS (see below). The seis- 
mic properties of these models arc very different from the 
observations. By comparing small separations between the 
oscillations, we confirm that both models are much older 
than the observed stars. 

The reason for such a great age for the system is that the 
observed luminosity of a Cen A is much greater than that of 
Q Cen B according to their masses. Because a Cen A evolves 
faster than a Cen B, an old age is required for a simultaneous 
agreement between the models and observations. However, 
we should also question the accuracy of the observed values 
(e.g. the radii; see below). 

Recently, del Peloso et al. (2006) determined the age of 
the galactic thin disc using Th/Eu cosmochronology from 
the analysis of selected stars, including the components of 
a Cen, as 8.8 Gyr. This age is very close to the age we 
found from the models of a Cen A and B for non-seismic 
constraints. 

According to seismic properties of its components, the 
a Cen system has a much smaller age than 8.88 Gyr; in fact, 
it is about 5.6-5.9 Gyr. With such an age, however, the lumi- 
nosities of Q Cen A and B can not simultaneously be fitted 
into the observed luminosities. Therefore, there is no con- 
straint on luminosity of a Cen B in the calibration of these 
models. Because of this strategy, the model luminosities of a 
Cen B based on the seismic properties are 15 percent larger 
than the observed luminosity. To remove this discrepancy, 
either the mass of a Cen B can be reduced to 0.914 Mq or 
its radius can be 7 per cent higher than the observed radius 
in KervcUa et al. (2003). However, the discrepancy may be 
a result of some processes, such as rotational mixing, which 
are not involved in the model computations. 

For all the models that also satisfy the seismic proper- 
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tics, in contrast to the models for non-seismic constraints, 
qb is greater than ua- If this is the case, then the mixing- 
length parameter is not a simple function of stellax mass 
(Yildiz et al. 2006) at all evolutionary phases, but is also 
a function of time. Therefore, we derive two expressions (/i 
and f2\ see Section 3.4) based on the structure of the con- 
vective zone of a given model and also compatible with the 
results in Yildiz et al. (2006) near the ZAMS. While one of 
the expressions gives qa as an increasing function of time, 
the other gives it as a decreasing function. Meanwhile, ac- 
cording to both expressions, as increEises. Such expressions 
may also explain why the mixing-length parameter varies 
for stars from phase to phase (see, for example, Ferraro et 
al. 2006) 

We also compute small separations between the oscil- 
lation frequencies from the alternative expressions given by 
Roxburgh & Vorontsov (2003) for models and the observed 
seismic data. In this case, we obtain more complicated situ- 
ation than in the classical seismic HR diagram. The small 
separation of models for a Cen A computed from the ex- 
pression given by Roxburgh & Vorontsov (2003) shows no 
explicit time dependence. 
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